Classical Analytical Mechanics and Entropy Production 
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The usual canonical Hamiltonian or Lagrangian formalism of classical mechanics applied to macro- 
scopic systems describes energy conserving adiabatic motion. If irreversible diabatic processes are 
to be included, then the law of increasing entropy must also be considered. The notion of entropy 
then enters into the general classical mechanical formalism. The resulting general formulation and 
its physical consequences are explored. 

PACS numbers: 45.20.-d, 45.20.Jj, 05.70.-a, 05.70.Ln 



I. INTRODUCTION 

In typical macroscopic classical mechanics problems, 
one deals with Hamiltonians describing a few degrees of 
freedom, say (P,Q) = (Pi, . . . , P„, Qi, . . . , Q n ). From 
a microscopic viewpoint, the Hamiltonian contains an 
enormous number of degrees of freedom, say N ~ 10 23 . 
To derive the macroscopic Hamiltonian from the mi- 
croscopic Hamiltonian requires the notion of thermody- 
namic entropy. In detail, suppose a microscopic Hamil- 
tonian H(P,Q,p,q) wherein (p,q) represent microscopic 
degrees of freedom. Discrete degrees of freedom are left 
implicit. For fixed classical macroscopic values for the 
reduced degrees of freedom (P, Q), one may consider the 
quantum microscopic energy eigenvalue problem [l[ 



H[P,Q, P 



d 



k= 1,2,. 



,0. 



(1) 



The thermodynamic entropy is determined by the micro- 
scopic degeneracy via the Boltzmann-Gibbs law Q 



S(E,P,Q) =k B In Sl(E,P,Q). 



(2) 



In principle, Eq.([2) may be solved for the energy in the 
form 



E = H(P,Q,S), 



(3) 



yielding the classical macroscopic Hamiltonian which is 
now an explicit function of entropy. Our purpose is to 
consider the physical consequences of including the en- 
tropy in classical canonical dynamics. 

In SecHIl we consider the adiabatic classical dynamics 
of macroscopic systems. Since the entropy S is uniform in 
time for adiabatic dynamics, the usual Hamiltonian and 
Lagrangian dynamics holds true, albeit with a macro- 
scopic Lagrangian L(Q,Q,S) which also explicitly de- 
pends on entropy. In adiabatic classical mechanics, both 
entropy and energy are conserved. 

In Sec lIIH we explore the consequences of diabatic pro- 
cesses wherein the entropy is increasing. In accordance 
with the first and second laws of thermodynamics the 
system energy is still conserved albeit some of the energy 
is converted into heat precisely defined by the increase in 



entropy. In detail, the Hamiltonian in Eq.([3]) along with 
friction force components (/i, . . . , f n ) determine the cou- 
pled equations of motion 



Pn 



Q a = 
dH(P, Q, S) 



dQ a 



dH(P, Q, S) 

dP a 
±fa(P,Q,S) , 



TS = -Q a f a > 0. 



(4) 



To assure an increasing entropy, the friction force com- 
ponents must be opposite to the velocity components. 
While the heating rate due to frictional forces TS is pos- 
itive, energy is nevertheless strictly conserved, i.e. 



E = H = 



(5) 



in virtue of Eqs.Q. 

In Sec lIVl well known cases of friction are discussed. 
The friction term is then generalized to a function of 
coordinates and velocities as f a — —R(Q,Q)v a , and the 
consequences are explored. 

In Sec|V]the example of the underdamped simple har- 
monic oscillator is considered. It is used as an pedaga- 
logical demonstration in the calculation purely thermod- 
unamic quantities such as temperature as fucntions of 
time for the formalism considered herein. 



II. ADIABATIC PROCESSES 

In an adiabatic process, the system entropy is uniform 
in time 



0. 



(6) 



In general, the first and second laws of thermodynamics 
may be written in the form Q 

dH = V a dP a - F a dQ a + TdS, 
dH{P, Q, S) 



V a (P.Q,S) ■■ 
F a (P.Q,S) = 
T(P,Q,S) 



dP a 

dHjP, Q, S) 

dQ a 
dH(P,Q,S) 
OS 



(7) 
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wherein the Einstein summation convention over index 
a is being employed, and T is the system temperature. 
Adiabatic macroscopic Hamiltonian dynamics is of the 
usual classical form 

Q a = V a {P.Q,S), P a = F a (P.Q,S) and S = 0. (8) 

To obtain the equivalent Lagrangian form of the equa- 
tions of motion, note that 



L= -H + V a P a . 
Eqs.®, © and © imply 

dL = P a dQ a + F a dQ a - TdS. 



(9) 



(10) 



The adiabatic equations of motion in Lagrangian form 
are 



d f dL 



dt \dQ°J QiS \dQ a jQ,s 



dL 



and S = 0. (11) 



The system temperature in the Lagrangian representa- 
tion follows from Eq. (fl0|) according to 



T = -(-) 



(12) 



Conservation of energy, 



\dQ a ; Q ,s 

E = 0, 



(13) 



follows directly from Eqs. lfTTj) . 

Both energy and entropy are conserved during a classi- 
cal adiabatic motion. When forces of friction are present, 
the energy is still conserved but the entropy increases 
with time in agreement with the first and second laws of 
thermodynamics. Let us see how this comes about. 



III. DIABATIC PROCESSES 

Our purpose is to derive the rule of conservation of 
energy in the presence of frictional forces. Starting from 
the Hamiltonian definition of energy E = H(P,Q, S), 
one may compute the time derivative 



(14) 



dP a dQ a ^ dS 
Employing Eq.Q for frictional forces yields 

E = Q a P a + (-P a + f a )Q a + TS = f a Q a + TS. (15) 
Energy conservation in the form 

E = 0, (16) 



is equivalent to the condition that frictional forces pro- 
duce heat according to the entropy rule 



TS=-f a Q a . 



(17) 



The components of the frictional force f a (a = 1 . . . n) 
will oppose the direction of motion Q a (a = 1 . . . n) 
giving rise to the second law of thermodynamics in the 
form 



S > 0. 



(18) 



The equations of motion in Lagrangian form with dia- 
batic processes can be found by extending Eq. lfTTj) to the 
case wherein S ^ 0. This may be done by directly adding 
the frictional forces to Lagrange's equations, which now 
read 



d f dL 



dL 



P a (Q a ,Q a ,S) 
F a (Q a ,Q a ,S) 



dt\d&) dQ a+fa ' 
dL(Q a ,Q a ,S) 



dQ a 
dL(Q a ,Q a ,S) 
dQ a 



More compactly, Eq (jl9|) asserts that 

Eq, Eq ~t~ fa ■ 



(19) 



(20) 



The energy expression from the entropy dependent La- 
grangian, 



E = E(Q,Q,S), 

E = Q«(^L)-L, 
\dQ a 



(21) 



may be shown, via Eqs. (fTU|) and (|19p . to be consistent 
with 



E = H(P, Q, S) 
dE = Q a dP a - F a dQ a + TdS. 



(22) 



Employing the energy conservation Eq. (|16[) and the di- 
abatic Lagrange's Eg . (f2"0|) , one again finds the heating 
rate 



TS=-f a Q a . 



(23) 



The above considerations show that the same rule for 
entropy production follows from both the Lagrangian and 
Hamiltonian formalisms. 



IV. THE GENERAL FRICTION FUNCTION 

We now examine in more detail the generalized fric- 
tional force f a - In theory and experiment, it is known 
that friction is generally a complicated phenomenon re- 
quiring regime-specific models and free phenomenology 
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parameters [4| [5[. For example, consider the differences 
between dry friction, viscous drag and atmospheric re- 
entry. Respectively, the dissipative force is constant, lin- 
ear and quadratic in velocity. Given this assortment of 
possibilities, a general approach must have enough flexi- 
bility to describe a wide range of phenomena while still 
conforming to certain physically intuitive notions. It is 
clear that any friction force would oppose the direction 
of velocity. It should also depend on the generalized co- 
ordinates such that the expressions involve tensors under 
arbitrary coordinate transformation; e.g. 

fa(Q,Q) = -Rab(Q,Q)Q b (24) 

satisfies these conditions. 

Examples of transport tensors include 

D a ab (Q) 

R ab = = (dry friction), 

\j9df(Q)Q d Q f 

Rab = Vab(Q) (viscous drag), 

j~f Pabcd{Q)Q C Q d i , \ /nt .\ 

Rab = , = (re - entry). (25) 

\j9fg{Q)Q f Q 9 

In all of the above cases, the entropy production from 
Eq. (fr7|) may be written 

TS = R ab (Q,Q)Q a Q b . (26) 

In the case of viscous drag, 

fa(Q-,Q) = -Vab(Q)Q b (27) 

which can be written as the derivative with respect to 
velocity of some function T(Q, Q), 

/.(M) = -!|*. (28) 
oQ a 

Hence the expression Eq. ((28|) can be written in the 
quadratic form 

T(Q,Q) = ^Vab(Q)Q a Q\ (29) 

which is otherwise commonly known as the Rayleigh dis- 
sipation function [|| Q [1]. 

V. THE UNDERDAMPENED OSCILLATOR 

Let us illustrate for a simple soluble model how the 
thermodynamic variables may change in time along with 
the purely mechanical variables. For this purpose we con- 
sider the underdamped simple harmonic oscillator. The 
Lagrangian is 

L = | (i 2 - c 2 * 2 ) , (30) 



with an entropy production of 

TS = 2m'yx 2 wherein 7 < uq. (31) 
The equation of motion for the underdamped oscillator, 
x + 2 7 i + loI = 0, (32) 
has the well known solution 



x — xoe 7 * cos(fM + (f)) wherein fl — y lu^ — j 2 . (33) 
If the harmonic oscillator has a constant specific heat 

C, 

mC = T^, (34) 

then Eqs. ([3"Tj) , (|33|) and (|34j) allow one to compute how 
the temperature of the oscillator will vary in time; i.e. 

T(t)=T +(jf) J*x(t>) 2 dt'. (35) 

For the case of zero initial phase <\> ~ 0, 

lim T(t) = T + = T + ^x 2 ,, (36) 

one finds that the temperature increase AT = T — To 
obeys 




C\J 
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FIG. 1: For an undamped oscillator with a quality factor 
Q = "f/ujo — 20.0, plotted as a "dotted" curve is the relative 
displacement x/xo as a function of time. Plotted as a solid 
curve is the temperature change AT/(T 00 — To) wherein T^, 
is the final equilibrium temperature after the oscillator stops 
moving. 
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where the oscillation quality factor Q = (ojq/j). We have 
plotted both the damped harmonic oscillator coordinate 
and the the temperature increase AT in FigQ]above as an 
illustration of how thermodynamic parameters as func- 
tions of time can be calculated similarly to computations 
of ordinary mechanical coordinates. 

VI. CONCLUSION 

The classical formalism of mechanics was extended to 
also include consideration of entropy. The results allowed 
for a revised form of Lagrange's and Hamilton's equa- 
tions which necessarily included energy dissipation due 
to frictional forces. It was shown that when the frictional 



forces vanished, that the entropy production was also 
and classical results were reproduced. However, in the di- 
abatic cases, the law of energy conservation gave rise to 
the conventional condition of entropy production in both 
the Lagrangian and Hamiltonian frameworks. A simple 
example was then considered and worked out in detail 
where the rate of entropy production was connection to 
the heating rate. The explicit connection between the 
motion of a body and the rise in temperature was shown 
in closed analytical form. A more complicated problem 
such as the entry into the atmosphere and burn of a me- 
teorite can in principle be treated by the methods here 
discussed although analytical solutions would appear un- 
likely. 
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